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Abstract. The neutron star low-mass X-ray binary EXO 0748—676 started a transition from outburst to quiescence 
in August 2008, after more than 24 years of continuous accretion. The return of the source to quiescence has 
been monitored extensively by several X-ray observatories. Here, we report on four XMM-Newton observations 
elapsing a period of more than 19 months and starting in November 2008. The X-ray spectra show a soft thermal 
component which we fit with a neutron star atmosphere model. Only in the first observation do we find a significant 
second component above ~3 keV accounting for ~7% of the total flux, which could indicate residual accretion. 
The thermal bolometric flux and the temperature of the neutron star crust decrease steadily by 40 and 10% 
respectively between the first and the fourth observation. At the time of the last observation, June 2010, we 
obtain a thermal bolometric luminosity of 5.6 x 10^'^ (d/7.1 kpc)^ erg s~^ and a temperature of the neutron star 
crust of 109 eV. The cooling curve is consistent with a relatively hot, medium-mass neutron star, cooling by 
standard mechanisms. From the spectral fits to a neutron star atmosphere model we infer limits for the mass 
and the radius of the neutron star. We find that in order to achieve self-consistency for the NS mass between 
the different methods, the value of the distance is constrained to be ^ 6 kpc. For this value of the distance, the 
derived mass and radius contours are consistent with a number of EoSs with nucleons and hyperons. 
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1. Introduction 

EXO 0748—676 is one of the most interesting transient 
low mass X-ray binaries (LMX Bs) observed to date . It 
was discovered with EXOSAT (jParmar et all Il986l) . It 
shows 8.3 minute X-ray eclipses every 3.8 2 hr, irregular 



dipping activity and type I X-ray bursts (jParmar et al, 
119861: iGottwaldet al.lflQS^ . 



Assuming Roche lobe overflow, a main- sequence com- 
panion , and a 1.4 Mq neutron star (NS) IParmar et al 



1986^ constrained first the inclination and companion 
mass to 75° < i < 83° and 0.45 Mq, respectively. From 



measurements of the radial velocity of the companion, 
using th e technique of Dopp l er to mography on optical 

spectra, Munoz-Darias et al.l ( 20091 ) derived a mass for n— — ■ 

the neutron star between 1 and 2.4 Mp. and a mass ra- (jCallowav et al.ll2010l) 

( 20091) 



NS mass of 1.27 Mq and constrained the mass ratio to 
0.075 < q < 0.105. For a 1.4 Mq NS this translates into 
a mass for the companion star of 0.11 < M2 < 0.15 Mq. 

Quasi-periodic oscillations (QPOs) were discovered 
during persistent emission at Hz and kHz freq uencies 
(jHoman et al.lll999l : iHoman fc van der Klidl2000h . Burst 
oscillations were first det ected at 45 Hz in an average o f 
38 type-I X-ray bursts bv lVillarreal fc Strohmaver (2004), 
which were interpr eted as the spin freque ncy of the neu- 
tron star. However, iGallowav et al. I (l2010^ reported oscil- 
lations in the rising phase of two type-I X-ray bursts at 
a frequency of 552 Hz and concluded that the latter is 
the spin frequency of the NS, while the 45 Hz oscillation 
may arise in the b oundary layer between the disc and the 
NS (lBalmanll2009l) or it could be a statistical fluctuation 



tio 0.11 < q < 0.28. Bassa et al 



used the same 
technique to analyse observations performed after a sig- 
nificant cessation of accretion activity in EXO 0748— 676 
was detected in 2008 and placed a lower limit on the 



From the analy sis of a photospher ic radius expansion 
(PRE) X-ray burst IWolfF et al.l (|2005l ) estimated the dis- 
tance of EXO 0748-676 to be 5.9-7.7 kpc Extending the 
analysis to several type-I X-ray bursts, IGallowav et al 
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(|2008a[) and lGallowav et all (|2008bl) estimated a distance 
of 7.4 ±0.9 kpc and 7.1 ±1.2 kpc, respectively, the lat- 
ter value taking into account the touch-down flux and the 
high-inclination of EXO 0748-676. 

EXO 0748—676 is the only source for which gravita- 
tionally redshifted absorption lines dur ing X-ray bursts 
have been reported ( Cottam et al. 200^. A gravitational 
redshift of z = 0.35 at the star surface was inferred from 
these measurements, which tra nslates to a mass-to- radius 
ratio of M/R = 0.152 M^g/km (jCottam et al.ll2002l ). This 



provides an empirical constraint on the equatio n of state 
(EoS) of dense, cold nuclear matter. Ozel ( 20061 ) used the 
value of the gravitational redshift, combined with an es- 
timate of the stellar radius, to infer a mass for the NS 
of M > 1.82Mq for a stella r radius of 12 km in the 
coordinate frame of the NS. Ranch et al . I (l2008l) regarded 
the line identifications from lCottam et a l. (2002) as highly 
uncertain from computations of LTE and non-LTE NS 
model atmospheres and derived instead a redshift of z 
= 0.24 with an alternative line identification and a NS 
radius of R = 12-15 km for the mass range M = 1.4- 
1.8 M0. The spectral features seen in the burst spectra 
from the initial data ( Cottam et al.l 120021) could not be 



reproduced in t he 68 burst spectra from a longer obser 



vation in 2003 (jCottam et al.ll2008l) . This could be due 



to chang es in the ionisation conditions in the NS pho- 
tosphere ( Cottam et al. 20081). However, it remains diffi- 



cult to reconcile the narrowness of the absorption lines 
detected in the initial data if they originate in the sur- 
face of the NS after the r ecent rapid spi n reported by 
Galloway etaP (|2010h fe.g. lLin et a"l]l2010l) . 

An alternative method to set constraints to the EoS of 
NSs consists in deriving their mass and radius from X-ray 
spectral fitting to the thermal emission component after 
accretion ceases. The theory behind this method is the 
so-called deep crustal heating of accreting NSs, according 
to which the nuclear ashes sink in the NS crust under the 
weight of newly accreting matter and with increasing pres- 
sure, undergo a sequence of nuclear transformations (pyc- 
nonuclear reactions) acc ompanied by heat deposition (e.g 
Haensel fc Zdunik 1990l ) . The heat gained via this process 
is lost during quiescent episodes, resulting in thermal emis- 
sion from the NS surface. The equilibrium temperature is 
set by the temperature of the NS core, which in turn de- 
pends on the long-term time-averaged mass accretion rate 
of the NS and the extent to which the core is able to cool 



via neutrino emission (jRutledee et al.ll2002l ). Thus, fitting 
the thermal component NS atmosphere models allows us 
to gain insight in the structure and composition of the 
NS crust and core, if a reliable distance estimate is avail- 
able. The significant cease of accretion of EXO 0748— 676 
in 2008 (see below) opened the possibility of constraining 
its mass and radius usi ng this method. 



Wolff et all (|2008al) reported on RXTE observations 
of EXO 0748-676 during August 2008 showing the low- 
est flux level, 6.8 x 10^^^ ergs cm^^ s^^ between 3 and 
12 keV, since the beginning of the RXTE mission in 1996. 
The flux level appeared near the detection threshold of the 



RXTE AU Sky Monitor (ASM) during September 2008. 
Subsequent Swift and RXTE observations of the source 
in October 2008 confirmed that accretion had largely 
ceased and EXO 748—676 was returning to quiescence 
Wolff et all (l2008bh . 

Following the significant cease of accretion, we 
proposed four long (30-100 ks) observations of 
EXO 0748-676 with XMM-Newton, three equally 
spaced during the first year and a fourth one in the 
second year. With the sensitivity provided by these ob- 
servations we aimed to obtain new constraints to the EoS 
of dense matter via 3 complementary methods: analysis 
of the heating and cooling curves and spectral fitting 
with NS atmosphere models. The cease of accretion of 
EXO 0748— 676 has been a l so mo i iitore d with Chandra 
and Swift (jPegenaar et al. l2009l 1201 1[ ). However, the 
shorter observations and the lower sensitivities of these 
observatories compared to XMM-Newton limit their 
capability to obtain significant constraints to the mass 
and radius of EXO 0748—676. In this paper, we report on 
our monitoring program with XMM-Newton (the first of 
these observations has been analysed also by IZhang et al. 
(120101) '). The set of four observations analysed here 
constitute the most sensitive dataset of EXO 0748— 676 
after its return to quiescent state. 



2. Observations 



The XMM-Newton Observatory djans en et alll200ll) in- 



cludes three 1500 cm X-ray telescopes each with an 
EPIC (0.1-15 keV) at the focus. Tw o of the EPIC imag- 
ing spectrometers use MOS CCDs ( Turner et al.l 12001 ) 
and one uses p n CCDs (IStriider et al.lj200ll) . The RGSs 
(0.35-2.5 keV, iDen Herder et al.l 120011 ) are located be- 
hind two of the telescopes. In addition, there is a co- 
aligne d 30 cm dia.meter Optical/UV Monitor telescope 
(OM, Im ason et al.l 120011) . providmg simultaneously cov- 
erage with the X-ray instruments. Data products were re- 
duced using the Science Analysis Software (SAS) version 
10.0. We present here the analysis of EPIC data, RGS 
data from both gratings and OM data. 

Table [1] is a summary of the XMM-Newton observa- 
tions reported in this paper. All the observations were 
performed using the EPIC Full Frame mode. In this mode 
all the CCDs are read out and the full field of view of 
the pn (MOS) is covered every 73 ms (2.6 s). The source 
count rate is well below the threshold of 6 (0.7) s~^ where 
pile-up effects become important for the pn (MOS) cam- 
eras. We extracted source events from a circle of radius 
between 37 and 47" centered on the PSF core and back- 
ground events from a circle with the same radius centered 
well away from the source. Ancillary response files were 
generated using the SAS task arf gen. Response matrices 
were generated using the SAS task rmfgen. Background 
subtracted light curves were generated with the SAS task 
epiclccorr, which corrects for a number of effects like 
vignetting, bad pixels, PSF variation and quantum effi- 
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Table 1. XMM-Newton observations of EXO 0748-676 since November 2008. T is the exposure time and C the EPIC 
pn 0.3-10 keV source persistent net count rate. In all cases the Full Frame Mode with the Thin filter was used for all 
the EPIC cameras, except for Obs 0560180701, for which the Medium filter was used. 



Observation 


Date 


Instrument 


T 


C 


ID 






(ks) 






^UUO INOVeillUcl U 


pn 


9Q 


U.U4:0 HZ u.uuu 










n 1 fin 4_ f) nnq 








9Q 


U.IDU zt u.uuz 


0605560401 


2009 March 17 


pn 


41 


0.525 ± 0.004 






MOSl 


43 


0.105 ±0.002 






MOS2 


42 


0.131 ±0.002 


0605560501 


2009 July 1 


pn 


100 


0.473 ± 0.004 






MOSl 


101 


0.117 ±0.002 






MOS2 


101 


0.115 ±0.002 


0651690101 


2010 June 17 


pn 


96 


0.457 ± 0.004 






MOSl 


98 


0.110 ±0.001 






MOS2 


98 


0.106 ±0.001 



ciency and accounts for time dependent corrections within 
a exposure, like dead time and GTIs. 

The SAS task rgsproc with the option spectrumbin- 
72m5=lambda was used to produce calibrated RGS event 
lists, spectra, and response matrices in an uniform wave- 
length grid space. We also chose the option keepcool^no to 
discard single columns that give signals a few percent be- 
low the values expected from their immediate neighbours. 
Such signals could be mis-interpreted as weak absorption 
features in spectra with high statistics. We generated RGS 
light curves with the SAS task rgslccorr. We used the 
SAS task rgscombine to add spectra from the same RGS 
and order of different observations and to add spectra from 
RGSl and RGS2 for one single observation in order to re- 
duce the statistical errors. 

The OM was operated in Image Mode with UVWl, 
UVM2 and UVW2 filters in Obs 1 and in Image+Fast 
Mode with U filter in Obs 2-4. In the Image mode the 
instrument produces images of the entire 17' x 17' FOV 
with a time resolution between 800 and 5000 s. In the 
Fast mode, event lists with a time resolution of 0.5 s from 
a selected 11" x 11" region are also stored. The SAS tasks 
omichain and omf chain were used to extract images and 
light curves of EXO 0748-676. 

3. Results 

3.1. X-ray lightcurves 

Fig. □ shows 0.3-10 keV EPIC pn EXO 0748-676 back- 
ground subtracted lightcurves with a binning of 50 s for 
each of the observations listed in Table [TJ Despite the low 
count rate of the source in these observations, the nar- 
row eclipses are clearly evident. Intervals used for spec- 
tral analysis (not affected by high background flares) are 
marked on the top of each panel. We inspected the light 
curves together with the hardness ratio (counts in the 



2-10 keV energy range divided by those between 0.3- 

2 keV) to search for variations in the out-of-eclipse emis- 
sion, which could reveal residual accretion. For instance, 
detection of dipping activity would indicate the existence 
of absorption in an ionised atmosphere and therefore the 
existence of an accretion disc. Detailed plots of the light 
curves are shown in Fig. [2l The eclipses ingresses and 
egresses are sharp. 

Outside of the eclipses the light curves do not show 
large deviations from the average count rate and the hard- 
ness ratio is consistent with being constant. As an exam- 
ple, the distribution of count rate between the flrst two 
eclipses in Obs 2 is consistent with a Poisson distribution 
of mean 0.54 s^^ (x^ of 21 for 13 degrees of freedom). 

3.2. X-ray spectra 

For each observation, we first extracted EPIC and RGS 
spectra excluding the times of eclipses. We rebinned all 
the EPIC spectra to over-sample the FWHM of the en- 
ergy resolution by a factor 3, and to have a minimum of 
25 counts per bin, to allow the use of the statistic. 
We used the RGS spectra with two different binnings: in 
Sections 13.2. 2ll3373l we rebinned the RGS spectra to over- 
sample the FWHM of the energy resolution by a factor 

3 and to have a minimum of 25 counts per bin, to be able 
to consistently use the statistic for both EPIC and 
RGS spectra. In Sect. 13.2.11 we rebinned the RGS spec- 
tra with different methods to be sensitive to narrow fea- 
tures (see below). For the EPIC (RGS) spectra we used 
the 0.3-10 keV (0.5-1.8 keV) energy band. Constant fac- 
tors, fixed to 1 for the EPIC pn spectrum, but allowed to 
vary for the EPIC MOS and RGS spectra, were included 
multiplicatively in order to account for cross-calibration 
uncertainties. 

W e performed spectral analysis using XSPEC ([Arnaudl 
1996h version 12.6.0. We used the photo-electric cross- 
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Fig. 1. 0.3-10 keV EPIC pn background subtracted lightcurves with a binning of 50 s. The narrow echpses are clearly 
evident. The thick horizontal lines mark the good time intervals used for spectral extraction. 



sections and abundances of Wilms et al. ( 2000f ) to account 
for absorption by neutral gas (the tbabs XSPEC model). 
Spectral uncertainties are given at 90% confidence (Ax^ = 
2.71 for one interesting parameter) and upper limits at 
95% confidence. 

3.2.1. RGS spectral analysis 

We examined the 0.5-1.8 keV (7-25 A) first and second 
order RGS spectra to constrain the iVn in the direction 
of the source and to search for the signature of narrow 
absorption and emission features. 

We first rebinned the RGS spectra for each observation 
to over-sample the FWHM of the energy resolution by a 
factor 3 and to have a minimum of 10 counts per bin. 

We could fit the RGS spectra of all the observations 
with a continuum consisting of a neutron star atmo- 
sphere model component, suitable for cooling NSs (see 
Sect. I3.2.2p . modified by photo-electric abso rption from 
neutral material. We obtained a C-statistic (ICashlll979l) 
of 283, 377, 612 and 561 for 242, 357, 591 and 530 degrees 
of freedom (d.o.f.) for Obs 1 to 4, respectively. The values 
of TVh in units of lO^i cm'^ were 1.6 ;'j';3, < 0.7, <0.5 and 



<1.5 for Obs 1-4, respectively. None of the observations 
showed narrow features. 

In order to increase the sensitivity to narrow lines 
and edges, we combined the RGS spectra of all the ob- 
servations and rebinned the combined spectra to over- 
sample the FWHM of the energy resolution by a fac- 
tor 3 and used the C-statistic. We could fit the spec- 
tra with the same model as the individual observations 
and obtained a C-statistic(d.o.f.) of 1651(1673) and TVh 
< 0.5 X 10^1 cm-2. 

To evaluate the effect of the absorption component 
in our fit, we substituted the tbabs component by the 
more recent tbnew component and re-fitted the com- 
bined spectra. We obtained a C-statistic of 1651(1673) 
and TVh < 0.4 X 10^1 cm-2. 

Finally, we rebinned the merged spectra imposing a 
minimum signal-to-noise ratio (SNR) of 5 per channel, in 
order to minimise the background effects and re-fitted the 
first order spectra (the second order spectra do not have 
enough SNR to be fitted). We obtained a C-statistic of 
115 (70) and A^h < 0.3 x lO^i cm^^. 

In summary, we obtained consistent values of A'^h for 
different models of interstellar absorption and rebinning 
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Fig. 2. Detail of 0.3-10 keV EPIC pn lightcurves with a binning of 50 s. 



methods. However, we obtained significantly different val- 
ues of TVfj when using different continuum components. 
This is expected due to the combination of limited en- 
ergy band of the RGS and low statistics of these observa- 
tions. As an example, substituting the neutron star atmo- 
sphere component by a power-law component in the fits to 
the merged spectra we obtained values of of 2.4 ± 0.8 
and < 1.6 X 10^1 cm-2 and C-statistic of 1647(1674) and 
103(71) for a rebinning with oversampling by a factor 3 
and with a minimum SNR of 5, respectively. Therefore, we 
conclude that the RGS spectra alone cannot constrain the 
value of Nu for these observations due to the poor statis- 
tics and the limited bandwidth, necessary to determine 
the continuum. 

We did not find any significant narrow feature in the 
RGS merged spectra. 

3.2.2. EPIC and RGS spectral analysis 

We next fitted the EPIC and RGS spectra simultaneously 
for each observation with a continuum consisting of a 
neutron star atmosphere component, to account for the 
thermal emission below ^3 keV, and a power-law com- 
ponent, to account for emission above ~3 keV, modified 



by photo-electric absorption from neutral material (model 
tbabs* (nsatmos+p o) in XSPEC). The nsatmos model 
(|Heinke et al.ll2006l) includes a range of surface gravities 
and effective temperatures and incorporates thermal elec- 
tron conduction and self-irradiation by photons from the 
compact object. It also assumes negligible (<10^ G) mag- 
netic fields and a pure hydrogen atmosphere. The param- 
eters of the nsatmos component are the NS mass and ra- 
dius, the effective temperature in the NS frame, the source 
distance and a normalisation factor, which parameterizes 
the fraction of the surface that is radiating. Throughout 
this work, we imposed a lo wer limit to t he rn ass of the 
NS of 1.27 Mq, as derived bv lBassa etall (|2009l) . Further, 
since NSs in LMXBs do not show evidence, such as spin- 
powered pulsations or cyclotron spectral features, for a 
strong magnetic field, we kept the normalisation of the 
neutron star atmosphere model fixed to 1, which corre- 
sponds to the entire NS surface emitting. 

We performed 3 fits per observation: first we left the 
mass of the NS as a free parameter during the fit and 
then we fixed its value to the canonical mass of 1.4 Mq 
and to 1.8 A/0 for comparison. We extrapolated the flux 
of the thermal component to the energy range of 0.01- 
100 keV to estimate the thermal bolometric fiux. The fits 
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Table 2. Best-fits to the 0.3-10 keV EPIC and 0.5-1.8 keV RGS spectra for each observation individually with the 
tbabs* (nsatmos+po) model, kpi is the normalisation of the power-law component in units of 10~^ ph. keV~ 



Ftot and Fpow are the 0.3-10 keV total and power law unabsorbed fluxes in units of 10 erg cm ^ s and 



Fhol 
NS 



the 0.01-100 keV nsatmos unabsorbed flux in the same units. The distance was fixed to 7.1 kpc and the nsatmos 
normalisation to 1 during the fits (see text). 



Observation 



[10^ 



[eV] 



M 

[M0I 



R 
[km] 



Fhol 
NS 



Xl (d.o.f.) 



0560180701 
0560180701 
0560180701 

0605560401 
0605560401 
0605560401 

0605560501 
0605560501 
0605560501 

0651690101 
0651690101 
0651690101 



0.70 ±0.11 
0.73 ±0.10 



0.70 



+0.13 



0.76 ±0.09 
0.78 ±0.08 
0.77 ±0.09 

0.48±0.07 
0.50±0.07 
0.50±0.07 

0.61 ±0.07 
0.64 ±0.07 
0.62 ±0.08 



118 ±4 
117±6 
1181^^^ 

112±3 
111±4 
111+^" 

-51 

114±3 
114±6 
114 t«g 

112±3 
110±5 
1 1 1 +^°° 



1.4 (f) 
1.8 (f) 
<2.4 

1.4 (f) 
1.8 (f) 
<2.4 

1.4 (f) 
1.8 (f) 
<2.1 

1.4 (f) 
1.8 (f) 
<2.2 



15.6±1.1 0.7±0.7 4: 



15.11?:^ 
15.7 1^:? 

15.6 ±0.9 

15.1 ±0.9 

15.2 tj j 
13.1 ±0.7 
12.1 ±1.0 
12.0 t^-J 

14.0 ±0.7 

13.3 ±0.9 
13.7 



0.7 



+0.1 
-0.8 



0.7±0.8 4; 



1.33 ±0.02 

1.34 ±0.02 
1.33 ±0.02 

1.01 ±0.01 

1.02 ±0.01 
1.02 ±0.01 

0.833±0.008 
0.839 ±0.008 
0.840±0.008 

0.831 ±0.008 
0.838 ±0.008 
0.833 ±0.008 



0.09 ±0.02 
0.09 ±0.02 
0.09 ±0.02 



1.36 ±0.01 

1.37 ±0.02 
1.36 ±0.02 

1.12 ±0.01 

1.13 ±0.01 
1.13 ±0.01 

0.915±0.009 
0.923 ±0.009 
0.923±0.009 

0.918 ±0.009 
0.928 ±0.009 
0.921 ±0.009 



0.85 (247) 
0.85 (247) 
0.86 (246) 

1.04 (288) 

1.03 (288) 

1.04 (287) 

1.09 (404) 
1.09 (404) 
1.09 (403) 

1.02 (466) 

1.02 (466) 

1.03 (465) 



were acceptable, with reduced xSi ^1 fo^ 250-450 
d.o.f., for all the observations. The inclusion of a power-law 
component improved the goodness of the fit only for Obs 1. 
The best-fit results are shown in Tabic [2] for a distance of 
7.1 kpc. The errors on the temperature and radius of the 
NS are significantly smaller when the mass of the NS is 
fixed. The value of A^h shows significant differences among 
observations. This clearly affects the accuracy with which 
we can measure the value of the NS radius and its effective 
temperature, and shows that a single observation is not 
sufficient to impose constraints on the mass and radius of 
the NS. 



3.2.3. Simultaneous fits for all the datasets 

In order to obtain the best possible constraints for the 
mass and radius of the NS with the XMM-Newton set of 
observations, we next fitted the EPIC and RGS spectra 
of all the observations simultaneously. We first assumed 
that the value of the absorption does not change among 
observations and tied the values of A^h, and mass and 
radius of the NS for all the observations. Since the contri- 
bution of the power-law component was only significant 
in the first observation (see Sect. I3.2.2p . we also tied the 
index of the power-law component for the simultaneous 
fit and allowed the normalisation of this component to 
vary from one observation to another. Similarly, we al- 
lowed the temperature of the NS atmosphere component 
to vary between observations to account for the expected 
cooling of the NS crust. Table [3] shows the results of the 
best-fit model for 3 values of the distance: 5.9, 7.1 and 
8.3 kpc. The fits are acceptable with a of 1.03 for 1419 
d.o.f. Fig. [3] shows the results of the fits for a distance 



of 7.1 kpc (note that we only show the EPIC pn spec- 
tra for clarity, but the fit included the EPIC MOSl and 
M0S2 and the RGSs spectra too). Fits of the thermal 
component to a blackbody model yielded similar xj, but 
implied an emitting area much smaller than a neutron 
star's surface, in agreement with the expectations (since 
the emergent spec tra at Tef f ^ 5 x 10 ^ K is very different 
from a blackbodv. iBrown et al. (1998)). Therefore, we do 
not discuss the blackbody fits further. Substituting the 
power-law component by a thermal Comptonisation com- 
ponent (model comptt in XSPEC) yielded a similar 
and the parameters of the NS atmosphere were not af- 
fected, i.e. the values obtained for the mass, radius and 
temperature of the NS are robust against changes of the 
model for the emission component above ^3 keV. 

Fig. m shows the contour plots obtained for the mass 
and radius of the NS for the three distances considered in 
Table [3] As we increase the distance of the source in the 
model, both the best-fit values of the mass and the radius 
and the allowed region for these parameters increase. 

Next, we investigated how the errors of the fit vary 
when some parameters are fixed or untied in the fit. First, 
since we have indications from the individual fits per- 
formed in Sect. 13.2.21 that the value of A^h may be dif- 
ferent among observations, we allowed this parameter to 
vary in the fits (Case 1 in Table |4|). Second, since the in- 
dex of the power-law component is not well constrained, 
we fitted the spectra fixing the index to the best-fit value 
(r = 0.24) and to 1 (Cases 2 and 3 in TableS]). Finally, we 
fixed the mass to the best-fit value of 1.78 Af© and to a 
more canonical value of 1.40 Af©, but allowed to vary the 
index of the power law (Cases 4 and 5 in Tabled]). 

The A^H adopts slightly different values when allowed 
to vary among observations. However, this small change of 
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Fig. 3. Left: EPIC pn spectra fitted with the model tbabs* (nsatmos+po) . Obs 1-4 are shown in black, red, blue and 
green, respectively. Right: Unfolded EPIC pn spectra. The colors are the same as in the left panel. 

Table 3. Best-fits to the 0.3-10 keV EPIC and 0.5-1.8 keV RGS spectra for all the observations simultaneously with 
the tbabs* (nsatmos+po) model (see caption of Table[2]for definitions). 



Observation A'^h 
[10^^ cm-^] 

0560180701 0.66 ±0.04 
0605560401 
0605560501 
0651690101 

0560180701 0.64 ±0.04 
0605560401 
0605560501 
0651690101 

0560180701 0.63 ±0.04 
0605560401 
0605560501 
0651690101 



kT 



[eV] 



118.2^2? 
111.9 

+24 
-20 
+24 
-20 



109.1: 

107.6: 



121.4 
114.9 



M 

[M©] 



R 
[km] 



Distance: 5.9 kpc 

1.52 ±0.5 11.8 1:^:2 0.24 ±0.7 1.715:2 

<0.27 
<0.24 
<0.43 

Distance: 7.1 kpc 



120.1^23 1-78 
113.6 t^f 
110.8 ±14 
109.5 tit 



+28 
-23 
+20 



112.0 lao 
110.7±18 



+0.4 
-0.6 



13.71:2.? 0.23 ±0.7 1.7 tli 
<0.26 
<0.23 
<0.42 



2.12 



HO. 4 
-0.8 



Distance: 8.3 kpc 
15.2 tJ:g 0.28 



1. 



, +3.0 
' -1.3 



<0.38 
<0.21 
<0.40 



Ftot 



1.33 ±0.01 
0.98 ±0.01 
0.885 ±0.009 
0.841 ±0.008 

1.32 ±0.01 
0.98 ±0.01 
0.880 ±0.009 
0.834 ±0.008 

1.31 ±0.01 
0.97 ±0.01 
0.876 ±0.009 
0.832 ±0.008 



0.09 ±0.02 

< 0.011 

< 0.007 

< 0.014 

0.09 ±0.02 

< 0.011 

< 0.007 

< 0.014 

0.08 ±0.02 

< 0.011 

< 0.007 

< 0.014 



1.36 ±0.01 
1.08 ±0.01 
0.98 ±0.01 
0.933 ±0.009 

1.35 ±0.01 
1.08 ±0.01 
0.97 ±0.01 
0.925 ±0.009 

1.34 ±0.01 
1.07±0.01 
0.97 ±0.01 
0.919 ±0.009 



xi (d.o.f.; 



1.03 (1419) 



1.03 (1419) 



1.03 (1419) 



A^H does not influence the values of the other parameters 
significantly. In contrast. Table 2] shows that the large un- 
certainty in the mass of the NS drives the uncertainties of 
the other parameters of the fit to unrealistically large val- 
ues. For example, we detect variations in the effective tem- 
perature of the NS atmosphere between observations as 
small as 1 eV, but the calculated uncertainty is more than 
an order of magnitude larger. The large uncertainty in the 
mass arises from the fact that there are three parameters 
required to specify a NS model fully: its m ass, radius and 
effective temperature ( Heinke et al. 20061) . Therefore, for 



a given effective temperature, there ar e several acceptable 



pairs Mjvs-Rats (see Appendix B of Heinke et al. 20061 



for details). The uncertainty in the index of the power- law 
component does not have a significant effect in the uncer- 
tainties of the other parameters. Therefore, in the studies 
of the variations of the temperature of the NS atmosphere 
among observations that follow we took into account the 
errors listed in Table |4] for Case 4. We also observe that 
changing the index of the power-law component from the 
best-fit value 0.24 to 1 has the effect of increasing signifi- 
cantly the normalisation of this component, although the 
contribution to the total flux remains unchanged. 

We investigated the decay shape of the thermal flux 
and the temperature of the NS atmosphere as inferred 
from the nsatmos parameters. We fitted the temperature 
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Fig. 4. Contour plots for the mass and radius of the NS for the fits shown in Table [3] for a distance of 5.9 (left), 7.1 
(centre) and 8.3 (right) kpc. The solid, dotted and dashed lines represent the 68, 90 and 99% probability contours. 
Examples of mass and radius predictions for representative equations of state wit h neutron and protons (blue line s) 
and including quarks, hyperons or kaons (red lines) are shown for comparison (see iLattimer fc Prakash|[200l l2007l ) . 



Table 4. Best-fits to the 0.3-10 keV EPIC and 0.5-1.8 keV RGS spectra for all the observations simultaneously with 
the tbabs* (nsatmos+po) model (see caption of Table [2] for definitions). 



Observation Nn kT^ff MR T fcpi Ftot Fpo„ xt (d-o.f.) 
[lO^'cm'^] [eV] [Me] [km] 

Case 1: 

0560180701 0.66to.o3 laO.Otj? l-77to;7 13.7^3.2 0.26 ±0.7 1.8t:ii 1.32 ±0.01 0.09 ±0.02 1.34 ±0.01 1.03(1416) 

0605560401 0.67 ±0.06 113.8124 <0-25 0.99 ±0.01 < 0.011 1.09 ±0.01 

0605560501 0.60 ±0.05 110.5125 <0-31 0.869 ±0.009 < 0.007 0.96 ±0.01 

0651690101 0.66 ±0.06 109.5 l^g <0.41 0.840 ±0.008 < 0.014 0.930 ±0.009 

Case 2: 

0560180701 0.64±0.04 120.1±16 l.rsto.l iS.Ttli 0.24(f) l-Tttl 1.32±0.01 0.09±0.02 1.35±0.01 1.03 (1420) 

0605560401 113.6 1^? <0.21 0.98 ±0.01 < 0.011 1.08 ±0.01 

0605560501 110.8 ±15 <0.14 0.880 ±0.009 < 0.007 0.97 ±0.01 

0651690101 109.3 ±21 <0.27 0.836 ±0.008 < 0.014 0.925 ±0.009 

Case 3: 

0560180701 0.65 ±0.04 119.6 ±16 1-79 to.l 13.8t'i';5 1(f) 5.3 ±1.3 1.31 ±0.01 0.08 ±0.02 1.34 ±0.01 1.04 (1420) 

0605560401 113.4 ±16 ' ' <0.55 0.98 ±0.01 < 0.009 1.08 ±0.01 

0605560501 110.6 <0.67 0.881 ±0.009 < 0.010 0.97 ±0.01 

0651690101 109.1 ±20 <0.85 0.837 ±0.008 < 0.013 0.927 ±0.009 

Case 4: 

0560180701 0.64±0.04 119.8±2.5 1.78(f) 13.7 ±0.5 0.24 ±0.7 l.TtH 1.32±0.01 0.09±0.02 1.35±0.02 1.03 (1420) 

0605560401 113.4 ±2.1 <0.26 0.98 ±0.01 < 0.011 1.08 ±0.01 

0605560501 110.5 ±2.0 <0.23 0.879 ±0.009 < 0.007 0.97 ±0.01 

0651690101 109.3 ±2.0 < 0.42 0.836 ±0.008 < 0.014 0.925 ±0.009 

Case 5: 

0560180701 0.62 ±0.04 121.3 ±1.7 1.40(f) 14.4 ±0.4 0.21 ±0.7 l.dtH 1.31 ±0.01 0.09 ±0.02 1.34 ±0.01 1.03 (1420) 

0605560401 114.8 ±1.6 <0.25 0.97 ±0.01 < 0.011 1.07 ±0.01 

0605560501 111.9 ±1.5 <0.21 0.874 ±0.009 < 0.007 0.96 ±0.01 

0651690101 110.7 ±1.5 < 0.40 0.831 ±0.008 < 0.014 0.92 ±0.009 



curve with two models: an exponential decay function of 
the form y{t) ~ ae^'*^*'^^/'^ ± 6, where a is a normalisa- 
tion constant, the start time of the cooling curve and 
T the e-folding time (timescale for the temperature to de- 
cay by a factor of e~^), and a power-la w component of 
the fo rm y(t) = A{t — Iq)^ . Following Deeenaar et al.l 
( 201ll) . we fixed to to 2008 September 5, which is be- 



tween the first non-detection by RXTE/PCA and the first 
Swift observation of the source. Taking into account only 
the XMM-Newton observations presented here, we found 
an e-folding time of 133.5 ±87.8 d, a normalisation of 
17.2 ±5.8 eV and a constant level of 109.1 ±2.2 for the 
exponential fit {xt of 0.06 for 1 d.o.f.). When fitting the 
decay curve with a power-law, we found parameter val- 
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Fig. 5. Decay of effective temperature (left panel) and thermal bolometric flux (right panel) with time for the XMM- 
Ncwton observations of EXO 0748— 676. The dashed and dotted lines show the fit to a exponential function and to a 
power law, respectively. The XMM-Newton, Chandra and Swift points are represented with black circles, dark grey 
squares and light grey diamonds, respectively. The XM M-Newton point s are t he values obtained for Case 4 in Table SI 
The Chandra and Swift points have been taken from iDeeenaar et al.l (j2011ll . 



ues of A = 141.0 ± 8.4 eV and B = -0.04 ± 0.01 {xl of 0.4 
for 2 d.o.f.). The results of these fits are shown in Fig. [5] 
We note that with only four points we cannot favour one 
of the above fits. However, the very small change of tem- 
perature and flux between Obs 3 and 4 (less than 2 and 
6% respectively) could point to a flattening of the cool- 
ing curve. If this is conflrnied with future observations, 
the current power-law fit could be discarded in favour of a 
broken power-law fit with a tentative break between days 
200 and 300. 

We next fitted the temperature curve taking into ac- 
count the XMM-Newton, Chandra and Swift values si- 
multaneously. We found an e-folding time of 220 ± 65 d, 
a normalisation of 14.0 ±1.4 eV and a constant level of 
107.6 ±1.5 eV for the exponential fit {xl of 0.39 for 19 
d.o.f.). When fitting the decay curve with a power-law, we 
found parameter values of A = 135.8 ±2.5 eV and B=- 
0.035 ±0.003 (xl of 0.51 for 20 d.o.f.). In summary, the 
best-fit values obtained when Chandra and Swift val- 
ues are taken into account are consistent within the errors 
with the values obtained only with the XMM-Newton val- 
ues. We note that we took t he Chandra and Swift values 
from Degenaar et al. ( 20 111) , who fixed the mass of the NS 
to 1.4 Mq and the distance to EXO 0748-676 to 7.4 kpc 
in their fits. Although we fixed the mass of the NS to its 
best-fit value of 1.78 M© and the distance to 7.1 kpc to 
obtain the XMM-Newton points. Table H] shows that the 
temperature values are consistent within the errors for a 
mass of 1.78 Mq (Case 4) and 1.4 Mq (Case 5). In addi- 
tion. Table [3] shows that increasing the distance from 7.1 
to 8.3 kpc causes significant changes in the values of the 
mass and radius of the NS but not in the effective tem- 
perature. This explains that we obtain consistent decay 
curves when the Chandra and Swift values are included 



in the fit. We attribute the small errors of the Chandra 
points compared to the XMM-Newton points (despite the 
poorer quality of the Chandra spectra) to the fact that 



iDegenaar et al.l ()2011f) fixed all the parameters of the fit 
except the effective temperature and the normalisation of 
the power-law component. Therefore, the added value of 
the Chandra and Swift points should be taken with cau- 
tion. 



3.3. OM data 

We extracted images from the UVWl, UVM2 and UVW2 
filter OM exposures for Obs 1. Based on the expected 
decreasing optical magnitude of the source, we chose to 
perform Obs 2-4 only with the U filter to maximise the 
sensitivity. OM data were not available for Obs 3 due to a 
technical error. For Obs 2 and 4 we extracted images and 
light curves for all OM exposures. 

The images show a source consistent with the posi- 
tion of EXO 0748-676 which fades from Obs 1 to 4. In 
Obs 1, we obtained an average magnitude of 18.4 ±0.1 
in the UVWl filter. The source was not detected in the 
UVM2 or UVW2 filters. For Obs 2 and 4 we obtained av- 
erage U optical magnitudes of 21.7 ±1.8 and 22.3 ±3.2, 
respectively. We note that EXO 0748— 676 is not detected 
in 4 out of 10 exposures in Obs 2 and in 14 out of 19 
exposures in Obs 4. 

Fig. [6] shows the detections in all observations as a 
function of phase. Clearly, the optical magnitude was still 
decreasing in Obs 4 compared to Obs 2. Unfortunately 
the errors are too large since the source is falling below 
detectability with OM. Therefore, we cannot determine 
the presence of a modulation with orbital phase. 
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Fig. 6. Detections of EXO 0748-676 with the OM moni- 
tor as a function of phase. 

4. Discussion 

We analysed four XMM-Newton observations of 
EXO 0748— 676 spanning 19 months and starting in 
November 2008, when a significant halt of accretion into 
the NS was detected. 

The light curves show a persistent level of emission 
without significant variations of intensity, except for the 
eclipse periods, which are clearly seen in all the observa- 
tions. The persistent count rate decreased steadily by 40% 
from the first to the fourth observation. 

All the observations show a soft spectral component 
between 0.3 and ~3 keV, which we fitted with a NS at- 
mosphere model. The bolometric flux of the thermal com- 
ponent decreased by 30% between the first and the last 
observation. In parallel, we observe a 10% decrease in 
the effective surface temperature of the NS from 120 to 
109 eV, which we interpret as the cooling of the NS crust 
towards thermal equilibrium with the core, after having 
been heated by accretion during the 24 year outburst. We 
detected only a 5% decrease in fiux and a 1.2% decrease 
in temperature between the last two observations, indicat- 
ing that the NS crust may be close to reach thermal equi- 
librium with the core. Fitting the inferred temperatures 
with an exponential decay plus a constant offset yields an 
e-folding time of 133 ± 88 days. 

In addition to the thermal component, the first obser- 
vation shows a component above ^2-3 keV that we fitted 
with a power-law component, and which represents ~ 7% 
of the total 0.3-10 keV fiux. The index of the power law is 
poorly constrained, with a value of 0.24±0.7. This value is 
much lower than the one found bv lPeeenaar et ah (2011) 
in a previous analysis of this observation, 1.7±0.5, but 
the fiux contribution to the total flux is similar. The in- 
dex found in Obs 1 is consistent within the errors with 
values fou nd for power-law com ponents in other cooling 
NSs (e.g. ICackett et aDl2010bl ). Although the origin of 
this component is still unknown, it has been suggested 



levels observed ( Campana et al. 19981) . The fact that we 



do not detect dips in the light curves indicates that resid- 
ual accretion onto the NS via an accretion disc is unlikely. 

We do not observe any significant contribution of the 
power-law component after 6 November 2008. We deter- 
mined upper limits to the contribution of this component 
to the total fiux of 1.1, 0.8 and 1.7% in the XMM-Newton 
observations on 17 March and 1 July 2009 arid 17 J une 
2010, respectively. In contrast. iDegenaar et al.l (|201ll ) re- 
port a significant changing contribution between 5 and 
15% of the power-law component to the total flux based on 
Chandra observations of the source in 10 February and 5 
June 2009 and in 15 April 2010 (note the proximity in time 
of the Chandra observations to the ones analysed in this 
work) . The sensitivity to detect and accurately model the 
power-law component is crucial to constrain the behaviour 
of the thermal component. In this sense, the higher effec- 
tive area of XMM-Newton compared to Chandra and the 
longer exposures make the observations presented in this 
work more suitable to determine the contribution of the 
power law to the total flux. However, the different power 
law fluxes could be explained if the source varied between 
the Chandra and XMM-Newton observations. 

We observe a significant decrease in the optical mag- 
nitude of the system. The most significant drop of in- 
tensity occurs after Obs 1 and is very small between 
Obs 2 and 4. This supports further the results from X-ray 
spectral fitting in the sense that any residual accretion 
prese nt in Obs 1 has disappeared in subsequent observa- 
tions. iHvnes fc JonesI ([2009) observed EXO 0748-676 be- 
tween November 2008 and January 2009 using Andicam 
and the SMARTS 1.3 m telescope. They found an aver- 
age magnitude of R = 22.4 and J — 21.3 for the optical 
counterpart and a periodicity consistent with the orbital 
period of the system, indicating that at the time of the ob- 
servations emission from the accretion disc and/or X-ray 
heated inner face of the companion star dominate the op- 
tical emission. We do not observe any significant modula- 
tion in Obs 2 and 4, indicating that the emission from the 
accretion disc has most likely disappeared and the inner 
face of the companion star has subsequently cooled down. 
However, given the large errors of the measurements in 
Obs 2 and 4, the existence of a small modulation due to a 
remaining temperature gradient between the illuminated 
and dark face of the companion cannot be ruled out. 

We attempted to constrain the value of the interstellar 
absorption with the high-resolution RGS spectra. This is 
important especially in the case that spectral variability 
in the power-law component is detected. Coupled varia- 
tions between the power-law index a nd the value of iVn 
were found for Aql X-1 and Cen X-4 (jCampana fc Stellal 



2003t I Campana et al.l 120041 ) which could be explained if 



the power-law arises as shocked e mission from the pulsa r 



wind and the infalling material (|Campana et al.l Il998f ) . 



Unfortunately, the absorption edges are not strong for 
these observations due to the low interstellar absorption 
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in the direction of the source, <10^^ cm~^, and therefore 
the value of iVn determined from the RGS spectra alone 
depends on the continuum used. However, since we only 
detect contribution of a power-law component in the first 
of the four observations, the variations in the value of A''h 
obtained are non-critical for our analysis. 

We do not find absorption lines in the high-resolution 
RGS spectra, in agreement with expectations. The detec- 
tion of lines from the NS surface in the spectra of qui- 
escent NSs wou ld allow the deterinination of the gravita - 



neutrino emission from the core. iRutledge et al.l (|2002[ ) 
simulated the thermal relaxation of the NS crust in the 
transient KS 1731—260 considering different impurity frac- 
tions and thus conductivities of the crust and the possible 
presence of " enhanced" core neutrino mechanisms such as 
direct Urea or pion condensation. Comparing the shapes 
of the curves in their Fig. 3 with Fig. [SJ we conclude that 
EXO 0748— 676 is consistent with having a NS crust with 
a high thermal co nductivity and low-im purity material. 



in agreement with Deeenaar et al. (l201ll) and consistent 



tional redshift ( Brown et all Il998l: iRutledee et all l2002l) with the other three cooling N Ss (IWiinands et all 12002 



and are therefore of great interest. However, given the 
high spin and inclination of EXO 0748— 67 6, the proba- 



bility of detecting su ch features is very low (jChang et al, 
120061: iLin et"alll2010h . 



4.1. Constraints from the cooling and treating curves 

The gradual decrease in thermal flux and neutron star 
temperature can be interpreted as the NS crust cool- 
ing down in quiescence after it has been heated dur- 
ing its long accretion outburst. We found an e- folding 
time for this cooling of 133 ± 88 d, a normalisation of 
17 ±6 eV and a constant level of 109 ±2 eV when fit- 
ting the temperature decay curve with an exponential 
law. The normalisation and the constant level are con- 
sistent with those inferred from a similar fit to Chandra 
(13.4 ± 0.2 eV and 107.9 ± 0.2) and Swift (17.2 ± 1.8 eV, 
106.2 ±2.5 eV) observations. In contrast, the e-folding 
time although consistent within the errors, is smaller 
than the one inferred from Cha ndra (192 ±10 d) an d 
Swift (257 ±100 d) observations (|Degenaar et al.ll201lh . 
When fitting the decay curve with a power-law, we 
find parameter values of A =141.0 ±8.4 eV and B=- 
0.04 ± 0.01, consistent with the values inferred from a simi- 
lar fit to Chandra (A = 134.4 ± 1 eV, B = -0.03 ± 0.01) and 
Swift (A =144.7 ±3. 8 eV, B =-0.05 ± 0.01) observations 
(jPegenaar et al.ll201l[) . 

The cooling curve of EXO 0748— 676 is much flatter 
compared to the currently cooling NSs KS 1731—260 and 
MXB 1658-298. KS 1731-260 shows cooling consistent 
with a power-law decay with ind ex -0.125 ±0.007 sti ll 
8 years after the end of outburst ( Cackett et al. 2010a[) . 
MXB 1658—298 has cooled instead following an exponen- 
tial decay with e-folding time of 465 ± 25 d, reaching a 
temperature of 54 ± 2 eV, which h as remained constant 
for ~1000 days ( Cackett et al.|[2008 ). In contrast, the cool- 
ing curve of the NS XTE J1701-462 can be fitted with an 
exponential law and sho ws an e-folding time of 120 



+30 , 
20 ^ 



(|Fridriksson et al.ll2010l ). very similar to the value found 



for EXO 0748-676 in this work. 

The decay time provides a meas ure of the thermal re 



laxation time of the NS c rust (e.g. IRutledge et al.l 12002 : 



Brown fc Cummind |2009|) . In particular, observing the 
drop in luminosity permits the measurement of the ther- 
mal timescale of the crust and the relative magnitude of 
the drop tells us about the presence or absence of enhanced 



2004l:lBrown fc Cummind l2009HRTdriksson et al.ll2010t ). 

From the cooling curve we infer that the NS crust 
is close to reach thermal equilibrium with the core. 
Therefore, we can compare our results with theoretical 
predictions of the quiesce nt thermal l u minos ities of ac- 
creting NS transients. Ya kovlev et al.l (j2004l) computed 



the quiescent bolometric thermal luminosity as a func- 
tion of long-term time-averaged mass accretion rate for 
several models of accreting NSs warmed by deep crustal 
heating. We calculated the average accretion rate dur- 
ing the last 11.5 years of outburst monitored with the 
RXTE/ASM. Considering only daily detections above 
3 (7, we obtained an average count rate of 1.20 s"^. We 
found an average ASM count rate of 0.68 s^^ during 
the XMM-Newton monitoring of EXO 0748-676 in 2003, 
for which a bolometric (0.1-100 keV) u nabsorbed flux of 



8.44 X 10 erg cm ^ s ^ was obtained (|Dfaz Trigo et al 



20061 : iBoirin et al. 2007t) . Using this flux as a conver- 
sion factor we obtain an average bolometric flux of 
1.49 X 10~^ erg cm~^ s~^ during the second half of the 
outburst. Comparing the thermal bolometric luminosity of 
5.6 X 10^^ (d/7.1 kpc)^ erg s"^ in June 2010 and the aver- 
age accretion rate during outburst w ith the heating curves 



calculated bv I Yakovlev et al.l (|2004[ ) (see their Fig. 3), we 
conclude that if the crust of EXO 0748— 676 has reached 
or is almost reaching thermal equilibrium with the core, 
a low to medium-mass NS with a near-standard cooling 
scenario is more likely than a high-mass star with signifi- 
cantly enhanced cooling. 

We note that the values of the thermal bolometric flux 
from which we calculate the thermal bolometric luminos- 
ity are similar to the values inferred from Swift data but 
higher tha n values from Chandr a data, as previously re- 
marked bv lPegenaar et aL ( 2011 ). This probably points to 
an offset in the flux calibration among these observatories. 

4.2. Constraints from spectral fitting to the thermal 
emission 

We attempted to constrain the mass and radius of the NS 
from spectral fit s to the obse r ved t hermal emission from 
the NS surface. 'Zhan g et al. I (I2OIOI) performed a similar 
analysis but only with the first of the observations pre- 
sented here. They found that fitting the thermal emis- 
sion with two different NS atmosphere models, nsgrav 
and nsatmos, yielded similar results. Therefore we only 



12 



Diaz Trigo et al.: XMM-Newton observations of EXO 0748— 676 in quiescence 



used the nsatmos model in our analysis. By fitting the 
four observations simultaneously we were able to reduce 
significantly the contours that define the ma ss and radius 
of the NS with respect to the analysis by IZhang et al, 
( 201(t 



Steiner et al. determined an empirical dense 



matter EoS from a heterogeneous dataset of six NSs. 
They found significant constraints on the mass-radius re- 
lation for NSs, and hence on the pressure-density relation 
of dense matter. For example, for NSs where the photo- 
spheric radius equals the NS radius, they found radii of 
11.0 and 10.6 km for NS masses of 1.5 and 1.8 Mq, re- 
spectively. For NSs where they allowed a photospheric ra- 
dius larger than the NS radius, they found radii of 11.8 and 
11.6 km for a NS mass of 1.5 and 1.8 Mq, respectively. We 
obtain a mass of 1.78 Mq and a radius of 13.7 km for a dis- 
tance of 7.1 kpc and of 1.53 Mq and 11.8 km for the lower 
limit of the distance of 5.9 kpc. Interestingly, our solution 
for a dis tance of 5.9 kpc ag rees very well with the EoS 
found bv lSteiner et al. ( 2010l ) and would imply a medium- 
mass NS, in agreement with the mass deri ved by compar- 
ison w ith the heating curves obtained by lYakovlev et al 



( 20041 ) (see Sect. 14. ip . In contrast, comparing the contours 
in Fig. 2] with predictions for the mass to radius relatio n 
for representative EoSs (e.g. iLattimer &: PrakashI 120011 ) . 
the solution for a distance of 7.1 kpc rules out most of the 
EoSs derived for an interior of nucleons and hyperons. 

There are several factors which can include uncertain- 
ties in the calculation of the distance for EXO 0748— 676. 
Firstly, EXO 0748— 676 is a high-inclination, dipping 
source, and as such there may be a changing partial ob- 
scuration of the NS and its expanding photosphere by 
the disc during the rad ius expansion of the type I bursts 
(|Gallowav et al.ll2008bl ). The ratio of the peak burst flux 
to the flux at touch-down (at the moment when the 
photosphere 'touches down' on to the NS) is expected 
to increase with increasing radius of the photosphere, 
with respect to low-inclination sources, and thus the dis- 
tance could be overest imated. Taking this into account. 
Galloway et al.l (|2008bl ) estimated a value for the distance 
to EXO 0748-676 of 7.1 ± 1.2 kpc. However, the value of 
the distance could be smaller if the obscuration at the 
touch-down flux were still underestimated. We note that 
in an ana lysis of all the dipping s ources observed by XMM- 



Newton ijPiaz Trigo et al.H2006|) , EXO 0748-676 showed 



the least ionised atmosphere and a relatively cool plasma 
was present at all phases, indicating that obscuration in 
this source was larger than in other dippers. This re- 
sult was confirmed by an analysis of 6 00 ks of high- 



resolu tion RGS spectra of EXO 0748 -676 (|van Peet et al 



200i). Secondly, the distance calculations from type I X 
ray bursts assume that at touch-down, the photosphere 



radius is equal to the NS radius. However, iSteiner et al. 



(|2010 l) argue that a photosphere with a radius larger than 
the NS provides internal consistency in their analysis and 
thus the assumption that the photospheric radius is equal 
to the NS radius is suspect. Further uncertainties arise 
from variations in the composition of the photosphere. 



the NS mass or variations in the typical maximum radius 
reached during the type I burst episodes (which affects the 
gravitationa l redshift, and hence th e observed Eddington 
luminosity) (jGallowav et al.ll2008bl ). 

Therefore, we conclude that major uncertainty in the 
current limits on the mass and radius of EXO 0748—676 
are driven by the uncertainty in the distance estimate. 
Moreover, in order to obtain consistency between the mass 
obtained from the heating curves and the spectral fits to 
the thermal component, the distance to EXO 0748— 676 
should be ^6 kpc. Alternatively, a low/medium mass NS, 
albeit with a rad ius of 14-15 km, is possible at ^7 kpc for 
harder EoS fe.g. iLattimer fc Prakashll2007l) . 

4.3. A revision on the distance estimate to 
EXO 0748-676 

All current distance estimates to EXO 0748—676 are based 
i n the analysis of thermonuclear (type-I) X-ray bursts 
Jjonker & Nelemans"2TOl IWolff et al'l 120051: lOzell 120061 : 
iGallowav et al. 2008a.b|). This method assumes that the 
peak flux for very bright bursts reaches the Eddington 
luminosity at the surface of the NS, at which point the 
outward radiation pressure equals or exceeds the gravi- 
tational force binding the outer layers of accreted mate- 
rial to the star. Once the Eddington flux is reached, the 
spectral evolution during the first seconds shows a maxi- 
mum in blackbody radius simultaneous with a minimum 
in color temperature, while the flux remains constant, in- 
dicating that the photosphere expands and the effective 
temperature decreases in order to maintain the luminos- 
ity at the Eddington limit. A large uncertainty in the 
theoretical Eddington luminosity arises from variations 
in the photospheric composition. A second large uncer- 
tainty arises from the fact that different bursts exhibit 
different peak fiuxes. These d ifferences may be cau sed 
by variable obscuration (e.g. iGallowav et al.l 120031) . If 
this is true, the uncertainty will be significantly larger 
for high-inclination sources, such as EXO 0748— 676, for 
which s trong absorption alon g the line of sight is common. 
Indeed , lAsai fc Dotani Hooi) found that the burst profiles 
of EXO 0748— 676 in the soft energy band are highly af- 
fected by the presence of a photo-ionised plasma, whose 
ionisation degree changes largely by the strong X-ray ir- 
radiation of the burst. Since absorption along the line of 
sight reduces the observed flux, bursts with higher flux 
will be more reliable as distance estimators since they are 
less affected by obscuration. Further uncertainties include 
the mass of the NS and the gravitational redshift. 

Up to now, three photospheric radius expansion bursts 
have been reported f rom RXTE on May 2004 an d 



June a nd August 2005 (| Wolff et al.l 120051: IGallowav et al 
2008al ). The burst on May 2004 shows the largest 
peak flux. Based on spectral analysis of this burst 



IWolff et al.l (l2005l) estimated a distance of 7.7 kpc 
(5.9 kpc) for a helium(hydrogen )-dominated burst pho- 
tosphere. Gallowav et al. ( 2008al ) used the three bursts 
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to infer a distance of 8.3/(l+X) kpc, where X is for 
solar abundances and 0.7 for hydrogen-rich materiaL 
Further, they constrained the distance to 7.4 kpc argu- 
ing that the bursts short rise times and low duration sug- 
gest ignition in a He- ri ch environment. We note that the 
Galloway et al.l (|2008al) obtain a peak flux for the bright- 
est burst of 4.7 x 10~^ ergs cm~^ s~ ^, 9% lower than the 
flux obtained by IWolff et al. ( 2005 ) due to the lower ef- 
fectiye area of LHEASO FT y5.3, used in the latter paper 
(|Galloway et al.ll2008al) . 

We estimate that Compton scattering in a highly 
ionised plasma (log ^ ^4) for a high-inclination source 
such as EXO 0748-676 can reduce the flux by 1 to 10% 
for a column density between 10^^ and lO^'^ cm~^. If the 
plasma has a lower degree of ionisation, the flux is reduced 
by a larger factor. Therefore, if the flux of the burst on 
May 2004 is underestimated by 1-10%, this translates into 
a distance of 6.8-6.5 kpc for a helium-dominated burst 
photosphere and of 5.2-5.0 kpc for a hydrogen-dominated 
burst photosphere. These yalues are consistent with an up- 
per limit to the distance of ^ 7 kpc, required by physical 
EoSs (see Fig.©. 

5. Conclusions 

We analysed four obseryations of EXO 0748—676 spanning 
19 months after a significant decrease of accretion on the 
NS was detected. We conclude that: 

- the emission aboye 3 keV decreased significantly 
within six months of the end of the accretion phase 

- the NS crust has a high thermal conductiyity and low 
impurity material 

- the NS crust is yery likely close to reach thermal 
equilibrium with the core 

- the NS has a low-medium mass and is cooling by 
standard mechanisms 

- a consistent solution between spectral fitting to the 
thermal spectra and the heating curves requires a distance 
of <,6 kpc to EXO 0748-676. 
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